Astronomy & Astrophysics manuscript no. ms 


©ESO 2012 


January 20, 2012 





Are OPERA neutrinos faster than light because of non-inertial 

reference frames? 

Claudio Germana 

INAF- Astronomical Observatory of Padova, Italy 
e-mail: claudio . germana@gmail . com 

Received / Accepted 

ABSTRACT 

Context. Recent results from the OPERA experiment reported a neutrino beam traveling faster than light. The challenging experiment 
measured the neutrino time of flight (TOF) over a baseline from the CERN to the Gran Sasso site, concluding that the neutrino beam 
arrives ~ 60 ns earlier than a light ray would do. Because the result, if confirmed, has an enormous impact on science, it might be 
worth double-checking the time definitions with respect to the non-inertial system in which the neutrino travel time was measured. 
An observer with a clock measuring the proper time r free of non-inertial effects is the one located at the Solar System Barycenter 
(SSB). 

Aims. Potential problems in the OPERA data analysis connected with the definition of the reference frame and time synchronization 
are emphasized. We aim to investigate the synchronization of non-inertial clocks on Earth by relating this time to the proper time of 
an inertial observer at SSB. 

Methods. The Tempo2 software was used to time-stamp events observed on the geoid with respect to the SSB inertial observer time. 
Results. Neutrino results from OPERA might carry the fingerprint of non-inertial effects because they are timed by terrestrial clocks. 
The CERN-Gran Sasso clock synchronization is accomplished by applying corrections that depend on special and general relativistic 
time dilation effects at the clocks, depending on the position of the clocks in the solar system gravitational well. As a consequence, TOF 
distributions are centered on values shorter by tens of nanoseconds than expected, integrating over a period from April to December, 
longer if otherwise. It is worth remarking that the OPERA runs have always been carried out from April/May to November. 
Conclusions. If the analysis by Tempo2 holds for the OPERA experiment, the excellent measurement by the OPERA collaboration 
will turn into a proof of the General Relativity theory in a weak field approximation. The analysis presented here is falsifiable because 
it predicts that performing the experiment from January to March/ April, the neutrino beam will be detected to arrive ~ 50 ns later than 
light. 
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1. Introduction 

! The OPERA experiment dThe OPERA collaboration! I201 ll) re- 
cently reported a neutrino beam traveling faster than light. The 
experiment measures the distributions of neutrino time emis- 
sion/detection over a baseline from the CERN to the Gran Sasso 
(CNGS) site. Data are collected within runs lasting for several 
montho In these data, the neutrino beam time of flight (TOF) 
turns out to be ~ 60 ns shorter than that calculated by taking 
the speed of light in vacuum. If confirmed, the result would have 
fundamental implications for modern physics. Basically, it re- 
views the fundamental post ulate by Albe rt Einstein: "Nothing 
can travel faster than light" (lEinsteinl 1905). 
Although the implications are intriguing, one should consider 
possible systematic errors that are as yet undisclosed. The 
OPERA team pointed out that the research for possible sources 
of systematic errors is still in progress. 

The scientific community has produced a plethora of works 
to interpret the puzzling neutrino speed result. These works were 
based on some sy stematic errors caused by the Earth's motion 
(lMonderenll20T 1|) or errors in synchronization of clocks by ei- 
ther a third clo ck moving into th e gravitational fi eld of Earth 
(IContaldil201ll) or GPS satellites (Ivan E lburg 201 1), which were 



later reviewed (iBesidal I20TTI Armando V.D.B. Assis 20110). 
However, several works were based on superluminal neutrino 
and its implications, therefore considering the OPER A neutrino 
speed result free of systematic erro rs (e.g iTamburini & Lavederf 
1201 It lLaveder & Tamburinill201 ll) . 

In this letter we focus on the question with respect to 
what a measure is made. The CNGS baseline has been mea- 
sured through GPS benchmarks to better than 20 cm in the 
European Terrestrial Reference Frame (ETRF; Colosimo et al. 
20110). However, from the point of view of a light ray trav- 
eling over CNGS (or a beam of particles), the TOF might 
need to be corrected for effects caused by the non-inertial ob- 
server ETRF, because the Earth rotates, orbits around the Sun 
and exhibits other motions (precession, nutation, polar motion). 
Nanosecond-precision measurements of OPERA are demanded 
to be in pulsar timing. To study the actual pulsar clock behavior, 
the time tagged by clocks running on Earth has to be corrected 
to that by an inertial observer. 

A reference frame that approximates an inertial one is the 
Barycentric Celestial Reference System (BCRS), which is non- 
rotating and located at the Solar System Barycenter (SSB). To 
turn the problem into the inertial frame BCRS, we use the 
Tempo2 software (Hob bs et all 120061: lEdwards et al.l I2006T) . a 



1 http://proj-cngs.web.cern.ch/proj-cngs/ 



2 http://vixra.0rg/pdf/l 1 10.0047v4.pdf 

3 http://operaweb.lngs.infn.it/Opera/publicnotes/notel32.pdf 
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software conceived to correct for non-inertial effects the time of 
arrival of photons coming from pulsars. This uses the most ac- 
curate planetary ephemerides to date (Standish 1998_3), in agree- 
ment with the IAU 2000 resolutions. 



2. From ETRF to BCRS: Tempo2 

Tempo2 is the software that best-models light rays traveling 
from a pulsar to the observatory, with an accuracy of ~ 1 ns. 
The software corrects the time of arrival of photons, as measured 
at the observatory, for general relativistic effects on both the 
photon-path and clocks running in a gravitational field, against 
which time is measured. In this framework Tempo2 is used to 
study how the inertial observer at SSB would see a signal travel- 
ing over CNGS. 

After transforming the ETRF coordinates of the sites into 
the International Terrestrial Reference System ITRS0, they are 
inserted into Tempo2. The software calculates the ephemerides 
of the vector components r(E, ?t/rc)> pointing from the SSB 
to the site S on the geoid at terrestrial UTC time turc- The 
software prints the components of the vector r'(T, turc) from 
the SSB to the geocenter F and those of the vector s(S, turc) 
from the geocenter to the site. The vector r'(T, tare) is con- 
structed in the barycentric frame BCRS, while s(S, t) refers to the 
Geocentric Celestial Reference System (GCRS), which is non- 
rotating and located at geocenter. However, s(S, turc) would 
have b een equivalent if it had been constructed in the frame 
BCRS (lEdwards et al.l2006l) . In s(L, t UT c) the software accounts 
for Earth rotation, precession-nutation, polar motion and irregu- 
larities in both polar motion and rotation. 
The vector r(2, turc) m BCRS reads 

r(X, t UTC ) = r x (S, t UTC )u x + r y (S, t UTC )u y + r-(__, t UTC )u z , (1) 

where nature) = ^(TJutc) + Si&turc) 0' = x,y,z) and u. v , 
Uy, u z are unit vectors in the barycentric frame BCRS. 

If the light ray is emitted at the CERN at the space-time coor- 
dinate r(CN, tuTc) an d is detected at Gran Sasso at r(GS, t' UTC ), 
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Fig. 1. Top: Timing residuals after subtracting the expected TOF 
from that of Eq. |@), neglecting the Einstein delay term. The sim- 
ulation is performed for the year 2010, expressed in Modified 
Julian Day (MJD=JD-2400000.5). The bin time of the temporal 
series is 10 minutes. Bottom: Zoom-in over five days. 



the two events are linked by the equation 



c(T'(GS,t' UTC )-T(CN,t UTC )) = \r(GS,t' UTC )-r(CN,t UT c)\ 

where c is the speed of light, t (Z, turc) is the Barycentric 
Coordinate Time (TCB) as a function of the site S and UTC 
time. The TCB is th e proper time e xperienced by the inertial ob- 
server at the SSB (Standish 1998). The TCB time is linked to 
UTC time by the relation 



The term A' E (GS, t' UTC ) - Ae(CN, turc) is the difference of the 
Einstein delay at the two clocks and at f UTC , turc - In the OPERA 
(2)experiment clocks were not corrected for A E , therefore they tick 
at UTC time, which is not a proper time. 



t(E, t UTC ) = turc + ls + A E (S, t UT c) , 



(3) 



where Is are le ap seconds to tie the UTC time to Terrestrial Time 
(TT) (~ 66 s; ISeidelmann & Fukushima1ll992|) . A E (L, t UTC ) is 
the Einstein delay, which accounts for relativistic effects on 
clocks running in a gravitational field. 
Expressing equation (f2| in UTC we obtain 



c (t' UTC - turc) = |r (GS, t' UTC ) - r (CN, t UTC )\ + 

- c (A' E (GS, t' UTC ) - A E (CN, t UTC )) , 



(4) 



from which the TOF (t' UTC - turc) of the signal in UTC proper 
time can be deduced. 

4 ftp://ssd.jpl.nasa.gov/pub/eph/planets/ioms/de405.iom.pdf 

5 For full accuracy it is suggested to use coordinates in ITRS. ETRF 
coordinates were transformed into ITRS through the online converter at 
http://www.epncb.oma.be/_dataproducts/coord_trans/index.php 



3. Results 

The geocentric coordinates of both the CERN and the Gran 
Sasso sites reported in Colosimo et al. (201 1) were inserted into 
Tempo2. These coordinates refer to the origin of the OPERA 
detector reference frame at Gran Sasso and to the target focal 
point at the CERN, 730534.61 m away. An additional baseline 
is that between the beam current tra nsformer (BCT) and the fo- 
cal po int at the CERN (743.391 m: iThe OPERA collaboration! 
1201 lb . Because the geocentric coordinates of the BCT were not 
found in the literature, we assume the CNGS baseline to be 
730534.61 m. An additional baseline of 743.391 m is equal to 
~ 20 arcsec on the Earth's surface and would induce a negligible 
timing correction in the relativistic term v • s/c 2 (see Sec. 13. U . 

3.1. Tying non-inertial clocks on Earth to the clock at BCRS 

We now study the TOF from equation © neglecting the term 
A' E (GS, t' UTC )-A E (CN, t UTC ), like in OPERA. Fig.[TJ(top) shows 
the TOF of the signal after subtracting the value one expects, 
730534.61 m/c=0.002436801 s. There is an excess of up to 
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Fig. 2. Einstein delay correction to tie the clock at CERN from 
UTC to TCB time after removing the fit first-order polynomial. 



~ 240 ns modulated by both an annual component and the 
Earth's rotation (Fig. Q] bottom). Fig. Q] shows that the neutrino 
TOF is never constant and depends on the epoch at which the 
measurement is performed. 

To understand what happens, we need to analyze the Einstein 
delay difference A' E (GS, f UTC ) - Ae(CN, tujc)- We print the 
Einstein delay for the clock at the CERN to tie UTC to TCB 
time (eq. 0). The amount of the correction at MJD=55 197.0 
(January 1, 2010) is up to ~ 16 s and grows with a slope of 
~ 1.5 x 10~ 8 s/s since MJD=43 144.0003725. The linear drift of 
~ 1.5 x 10~ 8 s/s takes into account the linear term of the Einstein 
time-dilation integral, and another term acc ounting for the grav- 
itatio nal redshift due to Earth's potential (llrwin & Fukushima 
Il999h . The Einstein integral accounts for the special relativistic 
time-dilation (v 2 /c 2 , v velocity of the geocenter) and the grav- 
itational redshift at geocenter (U/c 2 , U gravitational potential 
where the geocenter moves). Because v ~ 30 km/s, the special 
time-dilation is v 2 /c 2 ~ 10~ 8 s/s. 

Fig. |2] shows the behavior of the Einstein delay correction after 
fitting a first-order polynomial over five days and subtracting it0 
The overall linear behavior hides a diurnal modulation of some 
microseconds owing to Earth's rotation. The modulation origi- 
nates from the special relativistic time-dilation on the geoid with 
respect to the geocenter. The correction reads v • s/c 2 , with v 
the velocity vector of the geocenter and s the site position vector 
pointing from the geocenter to the site. Because |v| ~ 30 km/s 
and \s\ ~ 6300 km (radius of Earth), the corrections amount to 
up to ~ 2 yus modulated by Earth's rotation. 
The non-zero slope in Fig. [2] after removing the linear drift sig- 
nifies that the corrections hide other terms. The Einstein delay 
integral, refe rring to the motion of the geocen ter, has several pe- 
riodic terms (Seidelmann & Fukushima 1992). 

The inertial observer ties the UTC time of the clock at Gran 
Sasso to the TCB as well. The corrections to the clock at the 
CERN are subtracted from the Gran Sasso clock corrections to 
give the Einstein delay difference A' E (GS, t' UTC ) - Ae(CN, tyre)- 
Fig. [3] (top) shows the result. The behavior is equal to that seen 
in Fig. [TJ Hence, the cause for the residual on each single TOF 
measurement shown in Fig. [TJis that the UTC time is not cor- 
rected for relativistic effects through the Einstein delay differ- 
ence A' E - Ae, as equation (O would instead require. Correcting 
the UTC time for the Einstein delay implies null residuals in 



6 The coefficient of the polynomial obtained from the fit is 
~ 1.5 x 1CT 8 s/s, as expected. 
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Fig. 3. Top: Corrections that the non-inertial observer ETRF 
should apply to keep the clock at Gran Sasso synchronized with 
that at the CERN. Bottom: The difference behavior of the veloc- 
ity components parallel to the position vector s at the two sites 
for the year 2010. 



Fig-E as one would expect. 

Because the UTC time in OPERA is not corrected for - Ae, 
the two clocks at the CERN and Gran Sasso are not properly 
synchronized. To keep them synchronized, the non-inertial ob- 
server at ETRF has to apply the corrections in Fig. [3] (top) to the 
UTC time at Gran Sasso. 

Like Fig.Q](bottom), the corrections to keep clocks synchro- 
nized oscillate with a diurnal behavior. Consequently, events that 
are simultaneous in TCB time are simultaneous in UTC time as 
well only when the correction in Fig. [3] (top) is zero. This hap- 
pens whenever the special relativistic correction on the geoid 
with respect to the geocenter, v • s/c 2 , is equal at both sites. 
To verify that the behavior in Fig. [3] (top) is indeed caused by 
differences in the term v • s/c 2 between the sites, Fig. [3] (bottom) 
shows the difference of the dot product v ■ j between the CERN 
and Gran Sasso (j unit vector along s). The difference of the ve- 
locity component parallel to s is up to 3 km/s, modulated by both 
diurnal and annual components. Its shape reflects that of the rel- 
ativistic corrections in Fig. [3] (top). A difference in velocity of 
up to 3 km/s corresponds to a difference in timing corrections 
(S \v\ x |s|)/c 2 ~ 210 ns. The annual component comes from the 
modu lation over the seasons of the angle between v and s (see 
also lMonderenll201ll) . 
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3.2. Distribution of timing corrections 

To check whether the overall effect of the corrections in Fig. [3] 
completely vanishes, we can study their distribution. 
We divide the y-axis of Fig. [3] (top) into 10 ns long channels. 
For each channel we count the number of events, integrating 
over a typical OPERA run. Fig. [4] (top) shows the distribution 
over an OPERA run from April 29 to November 22, 2010 (from 
MJD=55314.0 to MJD=55522.0). It peaks at both the -240 ns 
and +160 ns channel, therefore it might induce an overall error 
on the clock synchronization of ~ -80 nfl 
To investigate farther out, we study the distribution for measures 
taken over an entire year, from January 1 to December 31, 2010. 
Fig. |4] (middle) shows that the distribution is symmetric with 
respect to channel zero, inducing a null effect on the clock syn- 
chronization. 

The difference between the two distributions in Fig. [4] (top- 
middle) might be explained in terms of the variation of the angle 
between v and s over the seasons (see also Fig. [3] bottom). The 
OPERA runs have always been carried out from April/May to 
November, i.e., they missed the period from January to March 
that would be needed to complete one angle cycle. Therefore, 
if the measures are taken only in a window of the entire Earth 
orbit, it might be possible that the clock synchronization is af- 
fected by this systematic error and the non-inertial observer on 
Earth would measure a TOF shorter by ~ 80 ns. Because he is 
unaware of non-inertial effects, the non-inertial observer should 
take care to carefully synchronize the clocks. Alternatively, the 
measurements could be integrated over one angle cycle: In this 
way, from the perihelion to the aphelion and back, the effect van- 
ishes. 

If we add the ~ +80 ns correction to the discrepancy claimed 
in the OPERA experiment (-60 ns), the result might change: We 
obtain a neutrino beam arriving later than light. 

4. Conclusions 

The analysis by Tempo2 suggests that clocks on Earth should be 
tied to the TCB time to ensure their synchronization. The TCB 
time is the proper time experienced by an inertial observer lo- 
cated at the solar system barycenter. The UTC time ensures syn- 
chronization only when the relativ istic correction v-s/ c 2 (dubbed 
differential special time-dilation; Edwa rds et al.l20 06) is equal at 
both clocks. Differences in this correction of up to ~ 240 ns are 
seen and are modulated by both the Earth's rotation and revolu- 
tion. 

The overall effect of a non-synchronization of clocks may 
vanish over one Earth orbit. Integrating the non-synchronized 
clock over months might induce an overall effect on the order of 
several tens of nanoseconds. 

If the analysis described in this letter holds for the 

OPE RA experiment, t he excellent measurement by 

iThe OPERA co llaboratio n! (1201 ll) wou l d turn into a proof 
of relativistic theories (Einstein 1905, 1916). The analysis 
might cast doubts on interpretations dealing with superluminal 
neutr ino that are, or are not, accounting for relativisti c effects 
(e.g. lTamburini & Lavederl201 it iLa veder & Tamburiniil l201 ll) . 

To check whether or not the present study holds, it might 
be interesting to see results from OPERA runs performed in a 
different period of the year than usual. Fig. [4] (bottom) shows 



7 The ~ -80 ns overall effect holds for t he last OPERA results release 
as well dThe OPE RA collaboration 201 1), because the run lasted from 
October 21 to November 7, 2011. 
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Fig. 4. Distribution of timing corrections. Top: Distribution cal- 
culated over a typical OPERA run (April to November 2010). 
It would imply an overall miss-synchronization of ~ -80 ns. 
Middle: Distribution over an entire year. Bottom: Distribution 
from January 1 to March 31. It would imply an overall 
miss-synchronization of ~ +50 ns. 



the distribution of the corrections for non-inertial effects from 
January-March. Again, it is not symmetric with respect to chan- 
nel zero, but shows an opposite result to that of the distribution 
on the top. From this result the neutrino beam should arrive ~ 50 
ns later than light. 
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